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Abstract 



The behavior of the electronic system of heavy fermion metals is considered. We 
show that there exist at least two main types of the behavior when the system is 
nearby a quantum critical point which can be identified as the fermion condensation 
quantum phase transition (FCQPT). We show that the first type is represented by 
the behavior of a highly correlated Fermi-liquid, while the second type is depicted by 
the behavior of a strongly correlated Fermi-liquid. If the system approaches FCQPT 
from the disordered phase, it can be viewed as a highly correlated Fermi-liquid which 
at low temperatures exhibits the behavior of Landau Fermi liquid (LFL). At higher 
temperatures T, it demonstrates the non-Fermi liquid (NFL) behavior which can 
be converted into the LFL behavior by the application of magnetic fields B. If the 
system has undergone FCQPT, it can be considered as a strongly correlated Fermi- 
liquid which demonstrates the NFL behavior even at low temperatures. It can be 
turned into LFL by applying magnetic fields B. We show that the effective mass 
M* diverges at the very point that the Neel temperature goes to zero. The B — T 
phase diagrams of both liquids are studied. We demonstrate that these B — T phase 
diagrams have a strong impact on the main properties of heavy-fermion metals such 
as the magnetoresistance, resistivity, specific heat, magnetization, volume thermal 
expansion, etc. 
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In heavy-fermion (HF) metals with strong electron correlations, quantum phase transitions at zero 
temperature may strongly influence the measurable quantities up to relatively high temperatures. 
These quantum phase transitions have recently attracted much attention because the behavior of 
HF metals is expected to follow universal patterns defined by the quantum mechanical nature of 
the fluctuations taking place at quantum critical points (see e.g. [1,2]). It is widely believed that 
the proximity of the electronic system of HF metal to quantum critical points may lead to non- 
Fermi liquid (NFL) behavior. The system can be driven to quantum critical points (QCPs) by 
tuning control parameters other then temperature, for example, by pressure, by magnetic field, or 
by doping. When a system is close to QCP we are dealing with the strong coupling limit where no 
absolutely reliable answer can be given on pure theoretical first principle grounds. Therefore, the 
only way to verify what type of quantum phase transition occurs is to consider experimental facts 
which describe the behavior of the system. Only recently, there appeared experimental facts which 
deliver experimental grounds to understand the nature of quantum phase transition producing the 
universal behavior of HF metals. 

It is the very nature of HF metals that suggests that their unusual properties are defined by a 
quantum phase transition related to the unlimited growth of the effective mass at its QCP. Moreover, 
a divergence to infinity of the effective electron mass was observed at a magnetic field-induced QCP 
[3-5]. We assume that such a quantum phase transition to be the fermion condensation quantum 
phase transition (FCQPT), an essential feature of which is the divergence of the effective mass M* 
at its QCP [6,7]. FCQPT takes place when the density x of system tends to the critical density 
xfc so that M* oc 1/r where r is a distance from the QCP, r = \x — xfc\- Such a behavior does 
not qualitatively depend on the system's dimensions and valid in both cases of two-dimensional 
(2D) and three-dimensional (3D) Fermi systems [8,9]. Beyond FCQPT the system possesses fermion 
condensation (FC) and represents a new state of electron liquid with FC [7,10]. As soon as FCQPT 
occurs, the system becomes divided into two quasiparticle subsystems: the first is characterized by 
quasiparticles with the effective mass M FC , while the second one is occupied by quasiparticles with 
mass M£. The quasiparticle dispersion law in systems with FC can be represented by two straight 
lines, characterized by the effective masses M FC and M£, and intersecting near the binding energy 
Eq. Properties of these new quasiparticles with M FC are closely related to the state of the system 
which is characterized by the temperature T, pressure P, or by the presence of the superconductivity. 
We may say that the quasiparticle system in the range occupied by FC becomes very "soft" and is 
to be considered as a strongly correlated liquid. Nonetheless, the basis of the Landau Fermi liquid 
theory [13] survives FCQPT: the low energy excitations of the strongly correlated liquid with FC arc 
quasiparticles, while this state can be considered as a quantum protectorate [6]. The only difference 
between the Landau Fermi-liquid and Fermi-liquid after FCQPT is that we have to expand the 
number of relevant low energy degrees of freedom by introducing a new type of quasiparticles with 
the effective mass M FC and the energy scale E [6]. It is possible to provide a consistent picture of 
high-T c metals as the strongly correlated Fermi-liquid [11]. 

When a Fermi system approaches FCQPT from the disordered phase, its low energy excitations 
are Landau quasiparticles which can be characterized by the effective mass M*. This mass strongly 
depends on the distance r, temperature and magnetic fields B [8]. At low temperatures, it becomes 
a Landau Fermi liquid with the effective mass M* provided that r > 0. This state of the system, 
with M* strongly depending on T, r and B, resembles the strongly correlated liquid. In contrast to 
the strongly correlated liquid, there is no energy scale Eq and the system under consideration is the 
Landau Fermi liquid at T — > 0. Therefore, this liquid can be called a highly correlated liquid. Such 
a highly correlated Fermi-liquid was observed in non-superconducting Lai.7Sr .3CuO 4 [8,12]. 
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In this Letter, we continue to show that within the framework of FCQPT it is possible to un- 
derstand the NFL behavior observed in different strongly and highly correlated Fermi liquids such 
as high-T c superconductors [11] and heavy-fermion metals. We apply the theory of fermion con- 
densation to describe the behavior of the electronic system of HF metals and to show that there 
exist at least two main types of the behavior. If the system approaches FCQPT from the disordered 
phase it can be viewed as the highly correlated electron liquid, and the effective mass M* depends 
on temperature, M* oc T -1 / 2 . Such a dependence of M* leads to the NFL behavior of the electronic 
system. The application of a magnetic field (B — B c ) > B*(T) oc T 3 / 4 restores Landau Fermi liquid 
(LFL) behavior. Here B c is a critical magnetic field. At (B — B c ) > B*(T), the effective mass 
depends on the magnetic field, M*(B) oc (B — B c )~ 2 / 3 , being approximately independent of the 
temperature at T < T*(B) oc {B - 5 C ) 4/3 - At T > T*(B), the T~ 1 / 2 dependence of the effective 
mass and the NFL behavior are re-established. At T — > 0, the system becomes LFL with the effective 
mass M* oc 1/r. When the system has undergone FCQPT, it becomes a strongly correlated electron 
liquid, and the effective mass behaves as M* oc 1/T leading to the NFL behavior even at low tem- 
peratures. The application of a magnetic field (B — B c ) > B*{T) oc T 2 restores the LFL behavior. 
At (B — B c ) > B*(T), the effective mass M*(B) oc (B — B^ x l 2 , being approximately independent 
of the temperature at T < T*(B) oc ^B - B c . At T > T*(B), both the 1/T dependence and the 
NFL behavior are re-established. We show that the effective mass M* diverges at the very point 
that the Neel temperature goes to zero. It is demonstrated that obtained B — T phase diagrams 
have a strong impact on the main properties of HF metals such as the magnetoresistance, resistivity, 
specific heat, magnetization, volume thermal expansion, etc. 

We start with the case of a highly correlated electron liquid when the system approaches FCQPT 
from the disordered phase. FCQPT manifests itself in the divergence of the quasiparticle effective 
mass M* as the density x tends to the critical density xfc, or the distance r — > [8,9] 

M* oc r oc -. (1) 

\x — Xpc\ r 

Since the effective mass M* is finite, the system exhibits the LFL behavior at low temperatures 
T ~ T*(x) oc \x — x FC \ 2 [8]. The quasiparticle distribution function n(p, T) is given by the equation 

Sn - £( p,T)-Mr)-Tlni^%^ = 0. (2) 



Sn(p,T) ' " ' n(p,T) 

The function n(p, T) depends on the momentum p and the temperature T. Here Q = E — TS — fiN 
is the thermodynamic potential, and /i is the chemical potential, while e(p, T), 

6E[n(p)] 



is the quasiparticle energy. This energy is a functional of n(p, T) just like the total energy E[n( 
entropy S'fn^)] and the other thermodynamic functions. The entropy S^nfjo)] is given by the familiar 
expression 

S[n(p)} = -2 J [n(p, T) In n(p, T) + (1 - n(p, T)) ln(l - n(p, T))] (4) 

which results from purely combinatorial considerations. Eq. (2) is usually presented as the Fermi- 
Dirac distribution 
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n(p,T) = <h + exp 



(e(p,T)-//) 
T 



(5) 



At T — > 0, one gets from Eqs. (2), (5) the standard solution n F (p,T — > 0) — > 0(pf — p), with 
e(p ~ p F ) — n — Pf(p — Pf)/M*, where pf is the Fermi momentum, 6{j>f — p) is the step function, 
and M* is the Landau effective mass [13] 

1 _ l de(p,T^0) 
M*~ p dp lp=PF ' [ ' 

It is implied that in the case of LFL M* is positive and finite at the Fermi momentum pp. As a 
result, the T-dependent corrections to M*, to the quasiparticle energy e(p), and other quantities, 
start with T 2 -terms being approximately temperature independent. The Landau equation relating 
the mass M of an electron to the effective mass of the quasiparticles is of the form [13] 

^ = ^ + /ft(p,P,.)V PI n(p 1 ) I |i 5 . (7) 

Applying Eq. (7) at T < T*(x), we obtain the common result 

M " = l-W)/3 - <8> 

Here Nq is the density of states of the free Fermi gas and Fl(x) is the p-wave component of the Landau 
interaction. At x — > x F c, the denominator in Eq. (8) tends to zero and one obtains Eq. (1). The 
temperature smoothing out the step function 9(p F — p) at pp = (x/37r 2 ) 1//3 induces the variation of 
the Fermi momentum Ap F ~ TM*/p F . We assume that the amplitude F L has a short range q <C p F 
of interaction in the momentum space. It is a common condition leading to the existence of FC and 
nearly-localized Fermi liquids [7,14]. If the radius is such that go ~ Ap F ~ TqM* /p F , corrections 
to the effective mass are proportional to T at T ~ T . Here T oc \x — x F c\ is a characteristic 
temperature at which the system's behavior is of the NFL type. On the other hand, at T*(x) T , 
we have q 3> T*(x)M* /p F) and the system behaves as LFL at T ~ T*(x), so that the corrections 
to the effective mass start with T 2 terms. We can also conclude that the transition region is rather 
large compared with T*(x) being proportional to T . 

In the case of T ~ To, we again can use Eq. (8) with F[(p F + Ap F ) ~ F F {p F ) + AAp F , where 
A oc dFl(x)/dx. Substituting this expansion of F F {p F + Ap F ) into Eq. (8) we find that 

M- ~ " « JL. (9) 
AAp F M*T v ; 

In deriving Eq. (9) we assumed that the system is close to FCQPT so that (1 — N Fl(p F )) <^ 
NqAAp f . We can say that at T ~ To, Ap F induced by T becomes larger then the distance r from 
FCQPT, Ap F > \p F c —p F \, where p F c corresponds to x F c- Solving Eq. (9) with respect to M*, we 
obtain [8] 

M*(T) oc -L (10) 



The behavior of the effective mass given by Eq. (10) can be verified by measuring the thermal 
expansion coefficient which is given by [15] 
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x f d(S/x ) 



3K \ dx 



Here, P is the pressure and V is the volume. Substituting Eq. (4) into Eq. (11), one obtains that in 
the LFL theory coefficient is of the order a(T) ~ MlT/p 2 F K. By employing Eq. (10), one obtains 
that at T ~ T [16] 

a(T) oc aVT + bT, (12) 

with a and b being constants. This result is in good agreement with experimental facts obtained in 
measurements on CeNi2Ge2 [17]. 

The application of magnetic fields B leads to Zeeman splitting of the Fermi level. As a result, two 
quasiparticle distribution functions with Fermi momenta p l F and p F appear, so that p F < p F < p 2 F 
and Ap F = (p 2 F — p F ) ~ p BM* jp F . Here p is the electron magnetic moment. In the same way 
Eq. (10) was derived, we can obtain the equation determining M*(B) [8]. The only difference is 
that there are no contributions coming from the terms proportional to Ap F , and we have to take 
into account terms proportional to (Ap F ) 2 . Assuming that the system is close to the critical point, 
we obtain 




M*(B) ~ M — . (13) 



At T ~ T*(x), Eq. (13) is valid as long as M*(B) < M*(x), otherwise we have to use Eq. (1). It 
follows from Eq. (13) that the application of magnetic fields reduces the effective mass. If there exists 
a magnetic order in the system which is suppressed by magnetic field B = B c0 , then the quantity 
(B — B c0 ) plays the role of zero field, and Eq. (13) has to be replaced by the equation, 

/ 1 x2/3 

m * ( Mb^) • (14) 

At high magnetic fields, we expect Eq. (14) to be invalid because Ap F becomes too large so that 
Ap F > go- In that case, the effective mass still depends on the magnetic field but the proportionality 
given by Eq. (14) is not preserved, and the dependence on the magnetic field becomes weaker. 

At elevated temperatures T ~ T , the effective mass starts to depend on both the temperature 
and the magnetic field. A cross over from the 5-dependent effective mass M*(B) to the T-dependent 
effective mass M*(T) takes place at a transition temperature T*(B) as soon as M*(B) ~ M*(T). 
This requirement and Eqs. (10) and (14) give that 

T*{B)o,{B-B cQ fl\ (15) 

At T > T*(x), Eq. (15) determines the line on the B — T phase diagram which separates the region 
of the LFL behavior taking place at T < T*(B) from the NFL behavior occurring at T > T*(B). 
At T < T*(B), the system behaves like LFL with the effective mass M*(B), and corrections to the 
effective mass start with T 2 terms. In accordance with the LFL theory, the specific heat c ~ 7T, 
with 

1 (B)ocM*(B)oc(B-B c0 )- 2 ^. (16) 
The resistivity p behaves as p = po + A(B)T 2 , where the coefficient 
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A(B) oc (M*(B)) 2 oc(B- 5 c0 )" 4/3 . (17) 

It follows from Eq. (16) and (17) that the Kadowaki- Woods ratio K = A/^ 2 [18] is conserved. All 
these results obtained from Eqs. (14-17) are in good agreement with experimental facts observed in 
measurements on the HF metal YbAgGe single crystal [19]. The critical behavior of the coefficient 
A(B) oc (B - B cQ f at B -> B c0 described by Eq. (14) with (5 = -4/3 is in accordance with 
experimental data obtained in measurements on CeCoIn 5 which displayed the critical behavior with 
(5 = -1.37 ±0.1 [4]. 

In the LFL theory, the magnetic susceptibility x M*/(l — Fq ). Note, that there is no ferro- 
magnetic instability in Fermi systems related to the growth of the effective mass, and the relevant 
Landau amplitude F a > —1 [14]. Therefore, at T < T*(B), the magnetic susceptibility turns out to 
be proportional to the effective mass 

X (B) oc M*(B) oc (B - B c0 )- 2 / 3 , (18) 

while the static magnetization M B (B) is given by 

M B (B)<xB X <x(B-B c0 y/ 3 . (19) 

At T > T*(B), as it follows from Eq. (10), Eq. (18) has to be rewritten as 

X (T) oc M*(T) oc -L (20) 

The behavior of x{B) and M B (B) as a function of magnetic field B given by Eqs. (18) and (19) 
and the behavior of x(T), see Eq. (20), are in accordance with facts observed in measurements on 
CeRu2Si2 with the critical field B c q — > [20]. 

Consider the system when r — > 0. Then its properties are determined by the magnetic fields B 
and the temperature T because there are no other parameters to describe the state of the system. 
At the transition temperatures T ~ T*(B), the effective mass depends on both T and B, while at 
T < T*(B), the system is LFL with the effective mass being given by Eq. (14), and at T > T*{B), 
the mass is defined by Eq. (10). Instead of solving Eq. (8), it is possible to construct a simple 
interpolation formula to describe the behavior of the effective mass over all the region, 

M*(B, T) oc F(B, T) = J (21) 

d(B - B c q) I 6 + c 2 f{y)VT 

Here, f(y) is a universal monotonic function of y — y/T/(B — B c0 ) 2 ^ 3 such that f(y ~ 1) = 1, and 
f(y <C 1) = 0. It is seen from Eq. (21) that the behavior of the effective mass can be represented 
by a universal function F M of only one variable y if the temperature is measured in the units of 
the transition temperature T*(B), see Eq. (15), and the effective mass is measured in the units of 
M*(B) given by Eq. (14), FM{y) = aF(B,T)/M*(B) with a being a constant. This representation 
describes the scaling behavior of the effective mass. As seen from Eqs. (19) and (21), the scaling 
behavior of the magnetization can be represented in the same way, provided the magnetization is 
normalized by the saturated value at each field given by Eq. (19) 

M *^ T ) o, \ (22) 



6 



where C3 is a constant. It is seen from Eq. (22), that magnetization is a monotonic function of y. 
Upon using the definition of susceptibility, x = dMs/dB, and differentiating both sides of Eq. (22) 
with respect to B, we arrive at the conclusion that the susceptibility also exhibits the scaling behavior 
and can be presented as a universal function of only one variable y, provided it is normalized by the 
saturated value at each field given by Eq. (18) 

X(B,T) ^ 1 n _ f(y) + ydf(y)/dy 

X(B) X l + c 3 f(y)y +ZC3y (l + c 3 f(y)yY " ^ 

It is of importance to note that the susceptibility is not a monotonic function of y because the 
derivative is the sum of two contributions. The second contribution on the right hand side of Eq. 
(23) makes the susceptibility have a maximum. The above behaviors of the magnetization and 
susceptibility are in accordance with the facts observed in measurements on CeRu 2 Si 2 [20]. Note, 
that the magnetic properties of CeRu 2 Si 2 do not show any indications of the magnetic ordering at the 
smallest temperatures and in the smallest applied magnetic fields [20], that is B c q — * in that case. 
As a result, we can conclude that the QCP is driven by the divergence of the effective mass rather 
then by magnetic fluctuations, and FCQPT is the main cause of the NFL behavior. We can also 
conclude that the Neel temperature is zero in this case, because the magnetic susceptibility diverges 
at T — > 0, as it is seen from Eq. (20). A more detailed analysis of this issue will be published 
elsewhere. 

Consider the case when the system has undergone FCQPT. Then, there exist special solutions of 
Eq. (2) associated with the so-called fermion condensation [7]. Being continuous and satisfying the 
inequality < n (p) < 1 within some region in p, such solutions no(p) admit a finite limit for the 
logarithm in Eq. (2) at T -> yielding [7] 

e(p)-/i = 0, if < no(p) < 1; p t <p <p f , (24) 

where e(p) is given by Eq. (3). At T = 0, Eq. (24) defines a new state of electron liquid with 
FC [7,10], which is characterized by a flat spectrum in the (pf — pi) region, and which can strongly 
influence measurable quantities up to temperatures T < Tf. In this state, the order parameter of the 

superconducting state «(p) = \J (1 — ^o(p)) n o(p) has finite values in the (pf — pi) region, whereas 
the superconducting gap A x — > in this region, provided that the pairing interaction tends to zero. 
Such a state can be considered as superconducting, with an infinitely small value of A l5 so that the 
entropy S(T = 0) of this state is equal to zero [6,7]. 

When pf — > pi — > pf the flat part vanishes, and Eq. (24) determines QCP at which the effective 
mass M* diverges and FCQPT takes place. When the density approaches QCP from the disordered 
phase, Eq. (24) possesses non-trivial solutions at x = xfc as soon as the effective inter-electron 
interaction as a function of the density, or the Landau amplitude, becomes sufficiently strong to 
determine the occupation numbers n(p) which delivers the minimum value to the energy E[n(p)}, 
while the kinetic energy can be considered as frustrated [6]. As a result, the occupation numbers 
n(p) become variational parameters and Eq. (24) has non-trivial solutions no(p), because the energy 
E[n(p)} can be lowered by alteration of the occupation numbers. Thus, within the region p^ < p < pj, 
the solution n (p) deviates from the Fermi step function np(p) in such a way that the energy e(p) 
stays constant, while outside this region n(p) coincides with n F (p) [7]. Note, that a formation of the 
flat part of the spectrum has been confirmed in Ref. [22-24]. 

At r > when the system is on the disordered side, that is «(p) = 0, and the density x moves 
away from QCP located at xfc, the Landau amplitude Fl(p = Pf,Pi = Pf,x) as a function of x 



7 



becomes smaller, the kinetic energy comes into a play and makes the flat part vanish. Obviously 
Eq. (24) has only the trivial solution e{p = pp) = H, and the quasiparticle occupation numbers are 
given by the step function, np(p) = 0(pf ~p)- 

At Ai — > 0, the critical temperature T c — > 0. We see that the ordered phase can exist only at 
T = 0, and the state of electron liquid with FC disappears at T > [6]. Therefore, FCQPT is not 
the endpoint of a line of finite-temperature phase transitions. This conclusion is in accordance with 
Eq. (2) which does not admit the existence of the flat part of spectrum at finite temperatures. As 
a result, the quantum to classical crossover upon approaching a finite-temperature phase transition 
does not exist. In the considered case, one can expect to observe such a cross-over at T ~ Tf. 
On the other hand, Ai becomes finite if we assume that the pairing interaction is finite, and the 
corresponding x — T phase diagram becomes richer. Moving along this line, we can consider the 
high-T c superconductivity as well, see e.g. [6,7]. 

At finite temperatures T <C Tf, the occupation numbers in the region {pf—pi) are still determined 
by Eq. (24), and the system becomes divided into two quasiparticle subsystems: the first subsystem 
is occupied by normal quasiparticles with the finite effective mass M£ independent of T at momenta 
p < pi, while the second subsystem in the (pf — pi) range is characterized by the quasiparticles with 
the effective mass M* FC {T) [6,25] 

M* FC ^ PF Pl^. (25) 

There is an energy scale Eq separating the slow dispersing low energy part, related to the effective 
mass M FC , from the faster dispersing relatively high energy part, defined by the effective mass M£. 
It follows from Eq. (25) that Eq is of the form [6] 

Eq ~ AT. (26) 

The described system can be viewed as a strongly correlated one, it has the second type of the 
behavior and demonstrates the NFL behavior even at low temperatures. By applying magnetic fields, 
the system can be driven to LFL with the effective mass [21] 

M*(B) oc 7 =L==. (27) 
V-d — n C Q 

In the same way as it was done above, we find from Eqs. (25) and (27) that a cross over from 
the -B-dependent effective mass M*(B) to the T-dependent effective mass M*(T) takes place at a 
transition temperature T*(B) 



T*{B)<x^{B-B c q). (28) 

Equation (28) determines the line in the B — T phase diagram which separates the region of the LFL 
behavior at T < T*(B) from the NFL behavior occurring at T > T*(B). The existence of the B — T 
phase diagram given by Eqs. (15) and (28) can be highlighted by calculating the resistivity and 
the magnetoresistance [8]. The resistivity, which at T > T*(B) demonstrates the NFL behavior, at 
T < T*(B), exhibits the LFL behavior, p = p + A(B)T 2 . The B-T dia gram of the dependence of 
the effective mass on the magnetic field can be highlighted by calculating the magnetoresistance. At 
(B — B c q) > B*(T), the magnetoresistance is negative, and at (B — B c q) < B*(T), it becomes positive. 
This behavior of both the magnetoresistance and the resistivity is in agreement with measurements 
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on YbRh 2 Si 2 [3], when the system exhibits the second type of the behavior, see Eq. (28), while 
CeCoIn 5 and YbAgGe demonstrate the first type of behavior consistent with that given by Eq. (15) 
[4,19]. 

At T < T*(B), the coefficients 7 oc M*(B), X (B) oc M*(B), and A{B) oc (M*(B)) 2 , and we find 
that the Kadowaki ratio K and the Sommerfeld- Wilson ratio R oc x(B)/ / y(B) are preserved due to 
Eq. (27). The obtained B — T phase diagram and the conservation of both the Kadowaki and the 
Sommerfeld- Wilson ratios are in full agreement with data obtained in measurements on YbRli2Si2 
and YbRli2(Sio.95Geo.o5)2 [3,5,17]. Taking into account Eqs. (11) and (25), we obtain that in the 
case of the two quasiparticle subsystems the thermal expansion coefficient a(T) oc a + bT + cVT, 
with a, b and c being constants. Here, the first term a is determined by the FC contribution, the 
second bT is given by normal quasiparticles with the effective mass M£, and the third cy/T comes 
from a specific contribution related to the spectrum £ c (p) which insures the connection between 
the dispersionless region (pf — occupied by FC and normal quasiparticles [16,25]. At finite 
temperatures, the contribution coming from the third term is expected to be relatively small because 
the spectrum £ c (p) occupies a relatively small area in the momentum space. Since at T — > 0, the 
main contribution to the specific heat c(T) comes from the spectrum £ c (p), the specific heat behaves 
as c(T) oc a\\/T + b\T, with a\ and 61 being constants. The second term b\T comes from the 
contribution given by FC and normal quasiparticles. Measurements for YbRh 2 (Sio.95Geo.o5)2 show a 
power low divergence of 7 = c/T oc T~ a with a = 1/3 [17]. This result is in a reasonable agreement 
with our calculations giving a = 0.5. At lower temperatures, the relative contribution of the first 
term aiVT becomes bigger, and we expect that the agreement will also become better. Now we 
obtain that the Griineisen ratio T(T) = a(T)/c(T) diverges as T(T) oc 1/VT [16]. This results is in 
good agreement with measurements on YbRh2(Sio.95Geo.os)2 [5,17]. 

As it follows from Eq. (27), the static magnetization behaves as Mb(B) oc \[B — B c q in accor- 
dance with measurements on YbRh 2 (Sio.95Geo.o5)2 [5]. We can also conclude that Eqs. (21), (22), 
and (23) determining the scaling behavior of the effective mass, static magnetization and the sus- 
ceptibility are also valid in the case of strongly correlated liquid, but the variable y is now given by 
y = T/y/B — B c0 , while the function f(y) can be dependent on (pf — Pi)/pF- This dependence comes 
from Eq. (25). As a result, we can obtain that at T < T*(B), the factor dp/dT oc A(B)T behaves 
as A(B)T oc T/(B — B c0 ), and at T > T*(B), it behaves as A(B)T oc 1/T. These observations are 
in good agreement with the data obtained in measurements on YbRh 2 (Sio.95Ge .o5)2 [5]. 

It is worthy to note that in zero magnetic fields the Neel temperature is zero at FCQPT, because, 
as it follows from Eq. (10), the effective mass tends to infinity at FCQPT and makes the susceptibility 
be divergent. On the other hand, if there is the magnetic order and the Neel temperature is not 
equal to zero, the effective mass is finite and there is no FCQPT. As soon as the magnetic order is 
suppressed at B — > B c0 , that is the Neel temperature tends to zero, the effective mass M*(B) — > 00 
as it follows from Eq. (14). If the system has undergone FCQPT, again at B — > B c0 , the Neel 
temperature goes to zero and M*(B) — > 00, see Eq. (27). If B c0 = 0, the effective mass diverges at 
T — > 0, see Eq. (25), and the susceptibility x tends to infinity being proportional to the effective 
mass. In this case, the Neel temperature is equal to zero as well. Therefore, one may say that the 
effective mass M* diverges at the very point that the Neel temperature goes to zero. 

A few remarks are in order at this point. To describe the behavior of heavy-fermion metals, we 
have introduced the system of quasiparticles, as it is done in the Landau theory of normal Fermi- 
liquids, where the existence of fermionic quasiparticles is a generic property of normal Fermi systems 
independent of microscopic details. As we have seen, at T Tf, these quasiparticles have universal 
properties which determine the universal behavior of heavy-fermion metals. One can use another 



9 



approach constructing the singular part of the free energy, introducing notions of the upper critical 
dimension, hyperscaling, etc., see e.g. [1,2]. Moving along this way, one may expect difficulties. For 
example, having the only singular part, one has to describe at least the two types of the behavior. 
We reserve a consideration of these items for future publications. 

In conclusion, we have shown that our simple model based on FCQPT explains the critical behav- 
ior observed in different HF metals. In the case of such HF metals as CeNi 2 Ge 2 , CeCoIn 5 , YbAgGe, 
CeRu2Si2, etc., the behavior can be explained by the proximity to FCQPT, where their electronic 
systems behave like highly correlated liquids. In the case of such HF metals as YbRh2(Sio.95Geo.o5)2 
and YbRli2Si2 the critical behavior is different. This can be explained by the presence of FC in 
the electronic systems of these metals, i.e., by the fact that the electronic systems have undergone 
FCQPT and behave as strongly correlated liquids. We have shown that the basis of the Landau 
Fermi liquid theory survives in the both cases: the low energy excitations of both strongly correlated 
Fermi-liquid with FC and the highly correlated Fermi-liquid are quasiparticles. It is also shown that 
the effective mass M* diverges at the very point that the Neel temperature goes to zero. The B — T 
phase diagrams of both the highly correlated liquid and the strongly correlated one have been stud- 
ied. We have shown that these B — T phase diagrams influence strongly the effective mass and such 
important properties of HF metals as magnetoresistance, resistivity, specific heat, magnetization, 
susceptibility, volume thermal expansion, etc. 
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my stay in Atlanta. This work was supported in part by the Russian Foundation for Basic Research, 
No 04-02-16136. 
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